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STABILITY  REGIMES  IN  A  HELICAL  QUADRUPOLE  FOCUSING 
ACCELERATOR-THEORY  AND  SIMULATION 


A  number  of  recent  high  current  spiral  line  or  recirculating  accelerator  configurations 
utilize  strong  focusing  fields.  These  fields,  consisting  of  a  helical  quadrupole  field  (or  stel- 
larator  field)  and  an  axial  guide  field,  increase  considerably  the  energy  mismatch  tolerance 
of  the  device  and  provide  confining  forces  against  the  beam  space  charge  forces.1,2  Two 
such  devices  are  the  modified  betatron  accelerator3  and  the  spiral  line  induction  accelerator 
(SLIA).4 

The  use  of  strong  focusing  fields  has  a  potential  difficulty  in  that  they  can  lead  to 
various  types  of  beam  instabilities.5-'  It  has  been  shown  that  the  electron  beam  centroid 
can  be  i)  orbit  unstable  independent  of  electromagnetic  waveguide  modes,  ii)  three- wave 
unstable  or  iii)  fully  stable.7  Here  we  present  the  dispersion  relation  for  the  instability  and 
analytic  conditions  for  each  of  the  stability  regimes  as  derived  in  Ref.  7.  Note  that  the 
stability  conditions  are  valid  only  in  the  limit  of  zero  beam  current.  Particle  simulations 
in  each  regime  will  be  presented  and  compared  to  both  the  analytic  stability  conditions 
and  to  numerical  solutions  of  the  dispersion  relation. 

DISPERSION  RELATION 

We  consider  an  electron  beam  progagating  in  an  external  magnetic  field  configuration 
consisting  of  an  axial  guide  field  B0,  and  a  helical  quadrupole  field  ( Bqx ,  Bqy), 

Bqx  =  -Bqkq(x  sin  kqz  -  y  cos  kqz),  Bqy  —  Bqkq{x  cos  kqz  -f  y  sin  kqz),  (la  —  6) 

where  Bqkq  is  the  quadrupole  gradient  and  kq  is  the  wavenumber  of  the  quadrupole  field. 

Here  we  quote  the  results  of  Ref.  7,  where  th-’  wave  equation  is  solved  simultaneously 
with  the  beam  dynamics  equation  to  obtain  a  yu-  rsion  relation.  In  the  analysis  it  was 
assumed  that:  1)  the  electron  beam  propagates  within  a  perfectly  conducting  cylindrical 
waveguide  of  radius  rg  (induced  image  charges  and  currents  were  included),  2)  the  beam 
radius  and  beam  centroid  displacement  are  small  in  comparison  to  the  waveguide  radius 
and  3)  the  electron  beam  is  monoenergetic  with  velocity  c0  and  relativistic  factor  ~)0  — 
(1  -  0l  )-1/2,  where  fi0  =  va/c. 

The  electromagnetic  fields  in  this  case  were  solved  for  in  terms  of  right-hand  circularly 
polarized  (RHCP)  and  left-hand  circularly  polarized  (LHCP)  waveguide  modes.  It  was 
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further  assumed  that  the  TE\\  mode  would  have  the  largest  growth  rate,  because  it  has 
an  electric  field  that  is  peaked  on  axis.  The  resulting  dispersion  relation  is 


WrW(WuWs  =  k2b  [(fc  +  kq  -  u>/v0f  D-  Wr  +  (k  -  uo/v0y  D+Wt  ]  , 


(2) 


where  kb  =  2(/&/17)^1 /(70(/xjjr2  -  7)J{(fuurg)),  Ib  is  the  beam  current  in  kA,  Mnr9 
is  the  smallest  positive  zero  of  Bessel  function  J[,  k  =  k+  =  k-  —  kq,  k±  are  the 
wavenumbers  associated  with  the  RHCP  and  LHCP  waves,  u>  is  the  radian  frequency, 
Wr  =  u>2/c2  —  k 2  —  and  W\  =  o;2/c2  —  (k  4-  kq)2  —  ^2j  are  contributions  from 
the  RHCP  and  LHCP  TE\\  waveguide  modes  respectively,  =  K2  —  (d2  +  d\)  ac¬ 
counts  for  the  two  stable  beam  modes,  Wu  =  K 2  —  (d\  —  d\)  admits  the  two  poten¬ 
tially  unstable  beam  modes,  d\  =  A'2  +  A'f /2,  d\  =  ((A'f  +  K{/2)2  -  (K*  -  A'|))'/2  , 
K  =  k  ~  u/v0  +  kqf 2,  A'i  =  K0  -  Jfe„  Kl  =  (K0  -  kq/2)kq/2  -  fc2,  K2  =  A’,fc„ 
A'0  =  |e|Z?0//3070m0c2  is  the  relativistic  cyclotron  wavenumber  associated  with  the  ax¬ 
ial  field,  Kq  =  \e\Bql j30~)0m0c2  is  the  relativistic  cyclotron  wavenumber  associated  with 
the  helical  quaclrupole  field,  k3  —  (2(Ij,/17)/(/3272r2  ))a/2  and  =  K2  ^  A'iA'  -  A’2. 
Equation  (2)  is  mathematically  identical  to  the  dispersion  relation  in  Ref.  6  with  the 
vertical  field  set  to  zero. 

STABILITY  REGIMES 

The  dispersion  relation,  Eq.  (2),  contains  i)  a  region  of  orbital  instability,  ii)  two  re¬ 
gions  of  three-wave  instability  and  iii)  two  regions  of  stability.  The  stability  conditions  were 
obtained  from  the  dispersion  relation  in  the  limit  of  zero  beam  current.1  They  delineate 
stable  and  unstable  regions  of  ( kq,K0 )  space  for  given  values  of  7,  rg  and  Bqkq. 

i)  Orbit  Unstable  Regime 

The  electron  beeam  is  both  orbit,  and  three- wave  unstable  when  (d j  -  d2)  0. 

which  gives  the  unstable  values  of  K0  : 


Ecrii,2  —  kq/2  2  A  q  <  A  0  ACri(,3  =  kqj  2  2I\q. 


(3) 
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ii)  Three- Wave  Unstable  Regimes 

Numerical  solutions  of  the  dispersion  relation  (2)  indicate  that  the  three-wave 
instability  occurs  when  the  RHCP  waveguide  mode  [or  the  LHCP  waveguide  mode] 
intersects,  in  the  (u?,k)  plane,  the  appropriate  beam  mode  given  by  Wu  =  0  and 
{d\  —  d\)  >  0.  For  Ka  <  Kcrit,2,  the  three- wave  interaction  is  unstable  (Region  I)  when 
the  RHCP  waveguide  mode  intersects  the  beam  line  u>/v0  =  (k  4-  kq/ 2)  4-  \/ d\  -  d\ 
^  k  4-  kq  (for  A'0  >  0)  or  ~  k  -f  kq  -  K„  (for  0  <  Ka  <  Kcrit (2).  For  K0  >  A'cri<,3, 
the  three-wave  interaction  is  unstable  (Region  II)  when  the  RHCP  waveguide  mode 
intersects  the  beam  line  oj/v0  =  (k  +  kqj  2)  —  \/  d\  —  d\  ~  k  4-  Ka.  The  two  three- wave 
unstable  regimes  are  those  portions  of  (kq,K0)  parameter  space  that  do  not  satisfy 
the  stability  conditions  below.  Note  that  while  the  three- wave  instability  is  present  in 
the  orbit-unstable  regime,  the  orbit  instability  dominates  in  that  regime. 

iii)  Three- Wave  Stable  Regime  for  K„  <  •> 

Stability  is  achieved  when  the  waveguide  cut-off  frequency  is  sufficiently 
large  that  intersection  with  either  of  the  beam  lines,  defined  by  Wu  =  0  cannot  be 
achieved.  For  Ka  <  KCrit,2i  the  beam  is  stable  if 

>  kq  4-  2  (dj  -  dl)1^  ,  (4) 

where  q  =  (4/(7„  —  2))1^2  .  For  kq  <  q^u,  and  defining  (  =  k*(  1  -f  8A'"//)  —  2/  and 
/  =  (kq  —  qv- u/2),  the  electron  beam  is  stable  for  the  following  situations: 

a)  for  /  >  0  and  C,  >  0,  the  stable  range  of  Ka  is  given  by 

Kcrit ,1  =  kq/ 2  -  C1/2/2  <  Ko  <  Kcrita,  (5a) 

b)  for  /  <  0  and  ^  >  0,  the  stable  values  of  K0  are 

K„  <  smaller  of  ( Kcrit,\ ,  Kcrit,2) ,  (56) 

c)  for  /  <  0  and  (]  <  0,  all  values  of  K0  <  KCrit,2  are  stable. 

i v )  Three- Wave  Stable  Regime  for  A„  >  Krrn  * 

The  three-wave  interaction  is  also  stable  when  the  RHCP  waveguide  mode  does 
not  intersect  the  beam  line  u>/r0  -  (k  +  kqj 2)  —  \/d\  —  d\  for  K0  Kcrit^.  This 

ry  ry  1  f  2 

occurs  when  qnn  >  kq  —  2  ( d\  —  d\)  ,  for  which  there  are  two  stable  cases: 


3 


a)  for  kq  <  q/J-u,  the  stable  range  of  K0  is 


h0  '>  hCrit,3i  (bo) 

b)  for  kq  >  qfi  11,  the  stable  range  of  Ka  is 

ho  >  hcrit,4  —  kq/2  -f  £*^/2.  (66) 

In  th**  limi*  of  "mall  quadrupole  gradient,  large  70  and  K0  >  KcriU3,  the  stability 

condition8  is  approximately  Ka  >  kq  —  Mii/7o- 

The  various  operating  regimes  are  illustrated  as  functions  of  kq  and  B0  in  Fig.  1, 
for  70  =  7,  rg  =  3  cm  and  quadrupole  gradient  =  200  G/cm.  Since  the  stability 

boundaries  are  obtained  in  the  limit  of  zero  beam  current,  the  area  of  the  two  stable 
regions  will  shrink  slightly  as  the  current  is  increased. 

It  is  clear  from  the  analysis  that  the  value  of  B0  >  0  required  for  stability  increases 
with  beam  energy.  To  operate  in  the  stable  regime  with  a  fixed  value  of  Ba  >  0,  both  the 
maximum  allowable  quadrupole  gradient  and  the  quadrupole  wavenumber  must  decrease 
for  increasing  beam  energy.  Figure  2  is  a  plot  of  maximum  quadrupole  gradient  versus 
70  for  a  fixed  guide  field  B0  =  5  kG,  rg  =  3  cm  and  two  different  values  of  quadrupole 
wavenumber:  kq  =  0.1  cm'1  ( — )  and  kq  =  0.05  cm-1(-  -  -). 

NUMERICAL  RESULTS 

For  the  present  numerical  study  we  use  the  ELBA9  code,  a  three-dimensional  particle 
code  which  simulates  a  beam  propagating  within  a  cylindrical  metallic  pipe.  The  full 
set  of  Maxwell’s  equations  along  with  the  full  relativistic  motion  of  the  beam  particles 
are  included.  The  initial  beam  parameters,  calculated  by  the  STELMAT10  code,  are 
matched  to  the  field  configuration  to  minimize  initial  oscillations.  The  ELBA  code  contains 
numerous  beam  diagnostics.  These  include  an  emittance  diagnostic11  that  was  developed 
for  beams  with  x-y  coupling. 

Growth  rates  were  measured  by  analyzing  the  TE\\  mode,  for  which  the  B-  and  Er 
components  may  be  “projected  out”  from  the  electromagnetic  spectrum  in  a  straightfor¬ 
ward  manner.  The  growth  rate  is  then  obtained  as  a  function  of  and  k  for  those  cases  in 


which  an  unstable  mode  grew  above  the  background  “noise”.  Because  our  simulation  takes 
place  in  a  coordinate  system  that  moves  with  the  beam  (r,  6,(  —  ci  —  r),  these  growth 
rates  are  T/c  —  Im(ujjc  —  k).  This  corresponds  to  the  theoretical  result  only  in  the  case 
that  Im(k)  —  0  as  was  assumed  in  numerically  solving  the  dispersion  relation  (2). 

Simulations  have  been  performed  for  parameters  typical  of  a  high  current  beam,  I/,  = 
10  kA,  70  =  7  and  normalized  RMS  emittance,  en,rm»  =  -158  cm  —  rad.  We  fixed  the 
waveguide  radius  at  rg  —  3  cm  and  the  quadrupole  gradient  at  Bqkq  =  200  G /cm.  In  a 
typical  run,  a  1  meter  beam  was  transported  over  10  meters  in  the  presence  of  the  external 
focusing  fields.  For  these  runs,  we  set  kq  =  0.5  cm~l  and  varied  Ba  so  as  to  sample  each 
of  the  unstable  regions  and  the  larger  of  the  two  stable  regions  of  the  stability  diagram 
(Fig.  1).  The  smaller  stable  region  with  B0  <  0  was  considered  by  setting  kq  =  0.09  cm-1. 
The  (kq,B0)  location  of  each  simulation  is  indicated  in  Fig.  1.  We  found  stable  behavior 
in  the  analytically  stable  regions  and  physically  distinct  behavior  in  each  of  the  unstable 
regions.  Numerical  results  (T/c,uj/c,k)  are  sumarized  in  Table  1  for  the  kq  —  0.5  cm-1 
runs.  Theoretical  values,  from  numerical  solutions  of  the  dispersion  relation,  are  also  given. 

Simulations  in  unstable  region  I  were  performed  with  Ba  —  0,  1  and  2  kG.  The  beam 
centroid  (x-component)  and  beam  envelope  (major  and  minor  radii  of  the  beam  cross 
section)  are  plotted  versus  time  for  a  fixed  position  within  the  beam  in  Fig.  3  for  the 
Ba  =  1  kG  case.  Here,  the  beam  develops  a  macroscopic  transverse  motion,  qualitatively 
similar  to  that  assumed  in  the  analytical  model  (pencil  beam,  rigid  displacements).  In 
fact,  the  wavenumber  of  the  beam  centroid  motion  for  each  of  these  runs  matched  the 
linear  theory  within  5%.  When  the  transverse  displacement  reaches  its  peak  (saturated) 
value,  the  cross-sectional  area  and  emittance  of  the  beam  increase  and  current  loss  sets  in. 
This  was  also  observed  in  the  Ba  =  2  kG  case,  where  saturation  occurs  at  a  displacement 
of  ~  0.1  cm  after  8  meters  of  propagation  and  is  followed  by  a  steady  loss  of  current  and 
rise  in  emittance.  Complete  saturation  was  not  observed  in  the  Bc  =  0  case  within  the 
length  of  the  simulation.  The  Fourier  spectrum  of  the  B-  component  of  the  TEi{  mode 
for  the  B0  =  1  kG  run  is  plotted  versus  k  in  Fig.  4  and  compares  favorably  to  the  growth 
rate  versus  k  plotted  in  Fig.  5. 

In  unstable  region  II,  we  found  saturation  of  the  instability  with  displacements  in  the 


5 


r 


1  to  2  mm  range  without  loss  of  current  or  growth  in  emittance.  Such  runs  were  performed 
at  Ba  =  4.5,  4.75  and  5.0  kG.  Plots  of  beam  centroid  (x-component)  and  beam  envelope 
are  shown  in  Fig.  6  for  the  Ba  —  5  kG  case.  The  Fourier  spectrum  of  the  B.  component 
of  the  TE ii  mode  at  ct  —  600  cm  is  plotted  in  Fig  7  and  compares  favorably  to  the  plot  of 
T/c  versus  k  from  the  dispersion  relation  (Fig.  8).  Figure  9  shows  the  growth  and  decay  of 
the  linear  component  of  the  mode  at  k  —  1.5  cm*1 .  Higher  wavenumber  modes  are  also  in 
evidence  but  are  within  the  range  of  unstable  wavenumbers  predicted  by  the  linear  theory. 
In  the  Ba  =  4.5  and  4.75  kG  cases,  higher  wavenumber  modes  appeared  that  were  outside 
the  linear  range  of  instability,  suggesting  that  higher  order  linear  or  nonlinear  modes  had 
been  excited.  Note  that  in  the  5  kG  case  the  dispersion  relation  gives  nonzero  growth  rates 
for  0.6  <  k  <  3.0  cm*1  and  for  5.7  <  k  <  9  cm*1.  Our  numerical  parameters  were  such 
that  wavenumbers  in  the  higher  range  were  not  properly  resolved. 

Simulations  of  orbit-unstable  configurations  were  performed  at  Ba  =  4.0  and  4.3  kG. 
These  points  in  parameter  space  reside  in  an  orbit-stable  region  of  the  stability  diagram, 
but  are  revealed  to  be  orbit  unstable  when  the  disperison  relation  is  solved  numerically. 
This  is  indicated  in  Fig.  10,  where  T  versus  k  is  plotted  for  B0  =  4.3  kG.  In  each  of  the 
two  simulations,  the  major  radius  of  the  beam  expanded  to  make  contact  with  the  wall 
within  less  than  1  meter  of  propagation,  resulting  in  a  severe  loss  of  beam  current  . 

Runs  were  also  performed  in  each  of  the  two  stable  regimes,  at  kq  =  0.5,  Ba  =  6.0  kG 
and  at  kq  =  0.09,  B0  =  —4  kG.  No  indication  of  instability  was  observed  in  either  case. 

DISCUSSION 

Two  items  among  the  results  above  deserve  comment.  These  are  1)  the  observed 
saturation  without  current  loss  in  three- wave  unstable  regime  II  and  2)  the  extension  of 
the  orbit  instability  into  tliree-wave  unstable  region  II  at  high  current. 

The  saturation  of  the  instability  in  three-wave  unstable  regime  II  without  emittance 
growth  or  current  loss  is  not  well-understood.  Unstable  region  II  differs  from  unstable 
region  I  in  that  the  instability  occurs  at  higher  wavenumbers  and  with  a  larger  bandwidth, 
which  may  allow  the  possiblity  of  mode  competition  as  the  instability  develops.  Also, 
at  short  wavelengths,  a  moderate  spread  in  the  parallel  energy  of  the  beam  may  disrupt 
the  instability.  This  saturation  mechanism  was  suggested  in  Ref.  8.  Because  the  group 
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velocity  of  the  TEn  mode  is  close  to  c  in  this  regime,  it  is  possible  that  an  electron 
trapping  mechanism  could  cause  saturation,  but  the  multiple  wavenumbers  observed  both 
in  the  radiation  field  and  in  the  beam  centroid  motion  make  this  unlikely.  It  was  not 
clear  from  either  the  analysis  or  the  simulations  which,  if  any,  of  these  phenomena  are 
responsible  for  the  observed  saturation. 

The  analysis  of  the  dispersion  relation  suggests  that  the  orbit  instability  is  centered  on 
a  region  of  parameter  space  where  the  cyclotron  wavenumber,  A'0,  is  approximately  equal 
to  kq/ 2  [see  Eq.  (3)].  For  such  values,  a  particle  gyrating  around  the  axis  experiences  a 
constant  radial  force  from  the  quadrupole  fields,  i.e.,  T(k  =  0)  ^  0.  At  high  current,  how¬ 
ever,  the  bandwidth  of  the  instability  in  three- wave  unstable  region  II  widens  to  encompass 
k  =  0  for  parameters  which  lie  outside  the  orbit-unstable  regime  in  the  stability  diagram 
of  Fig.  1  (recall  that  the  stability  diagram  is  correct  only  as  /(,  — »  0).  This  effect  can 
occur  in  three-wave  unstable  region  I,  but  is  less  significant  because  of  the  narrow  unstable 
bandwidth  in  that  region.  In  fact,  the  unstable  region  expands  from  2.2  <  Ba  <  3.8  IcG, 
as  suggested  by  the  stability  diagram,  to  2.2  <  Ba  <  4.3  kG  (see  Fig.  10).  This  exten¬ 
sion  of  the  orbit  unstable  regime  into  three-wave  unstable  region  II  is  the  only  significant 
departure  from  the  stability  diagram  at  high  current. 

CONCLUSIONS 

In  this  study,  we  have  used  electromagnetic  particle  simulation  to  examine  the  physics 
of  the  three-wave  instability.  We  first  presented  the  dispersion  relation  of  Ref.  7  for  trans¬ 
verse  oscillations  of  an  electron  beam  propagating  within  a  drift  tube  in  the  presence  of 
external  solenoidal  and  helical  quadrupole  focusing  fields.  We  also  presented  analytic  sta¬ 
bility  conditions  from  this  dispersion  relation  in  the  limit  of  zero  beam  current  and  have 
found  that  these  stability  conditions  are  good  predictors  of  particle  simulation  results,  even 
at  currents  as  high  as  10  kA.  In  accordance  with  the  theory,  we  find  that  as  we  vary  the 
external  field  parameters,  the  simulations  show  two  stable  regimes  of  beam  propagation, 
two  three-wave  unstable  regimes  and  an  orbit-unstable  regime.  The  only  significant  de¬ 
parture  from  the  stability  conditions  at  these  high  currents  had  the  effect  of  extending  the 
orbit  instability  into  three- wave  unstable  regime  II.  Numerical  solutions  of  the  dispersion 
relation  were  found  to  be  in  qualitative  agreement  with  particle  simulation  results  in  all 
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cases,  although  the  simulations  gave  lower  growth  rates  in  many  cases.  Additionally,  we 
have  observed  saturation  of  the  instability  in  both  three-wave  unstable  regimes.  This  was 
associated  with  emittance  growth  and  the  onset  of  beam  loss  in  three-wave  unstable  region 
I,  but  was  without  these  effects  in  unstable  region  II.  The  saturation  mechanism  in  the 
latter  case  is  not  clear,  but  may  be  related  to  the  short  wavelengths,  wide  bandwidths  and 
possible  mode  competition  which  occur  in  three-wave  unstable  region  II,  but  are  absent  in 
region  I. 
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Table  1 


THEORY  SIMULATION 


B0(kG) 

r/c 

u>/c 

k 

r/c 

IjO /c 

k 

0 

0.031 

0.71 

0.24 

0.005 

0.64 

0.17 

1.0 

0.033 

0.78 

0.38 

0.015 

0.76 

0.37 

2.0 

0.043 

1.03 

0.72 

0.033 

0.81 

0.49 

4.0 

orbit-unstable 

orbit-unstable 

4.3 

orbit-unstable 

orbit-unstable 

4.5 

0.023 

0.98 

0.82 

0.025 

1.06 

0.81 

4.75 

0.016 

1.29 

1.16 

0.013 

1.06 

1.09 

5.0 

0.012 

1.74 

1.65 

0.011 

1.87 

1.50 

6.0 

stable 

stable 

Theoretical  values  of  (u >/c,k)  for  modes  with  peak  linear  growth  rates,  r/c,  and  cor¬ 
responding  numerical  values  for  these  unstable  modes.  Growth  rates,  frequencies  and 
wavenumbers  are  in  units  of  cm-1.  Note  that  for  cases  where  multiple  unstable  modes 
were  observed  ( Ba  =  4.5,  4.75  and  5  kG),  only  the  mode  corresponding  to  the  linear 
growth  rate  was  reported. 
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Fig.  1.  Plot  of  the  various  operating  regimes  for  To  =  7,  rg  =  3  cm  and  quadrupole  gradient 
Bqkq  =  200  G/cm. 


0  2  4  6  8  10  12  14 

ct  (m) 


Fig.  3.  The  major  radius  of  the  beam  cross  section  (solid  line),  the  minor  radius  (dashed 
line)  and  the  beam  centroid  x-component  (dotted  line)  are  plotted  versus  time  for  the 
Ba  =  1  kG  case. 
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( ar  b.  units) 


i.o 


Fig.  4.  The  Fourier  spectrum  of  the  Bz 
the  B0  —  1  kG  case. 


f/c  (cm 


0.030 


k  (cm  ') 

Fig.  5.  The  temporal  growth  rate,  as  calculated  from  the  dispersion  relation,  is  plotted  versus 
k  for  the  B0  —  1  kG  case. 
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ct  (m) 

Fig.  6.  The  major  radius  of  the  beam  cross  section  (solid  line),  the  minor  radius  (dashed 
line)  and  the  beam  centroid  x-component  (dotted  line)  are  plotted  versus  time  for  the 
Ba  =  5  kG  case. 
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Fig.  7.  The  Fourier  spectrum  of  the  Bz  component  of  the  TEU  mode  is  plotted  versus  k  for 
the  B„  =  5  kG  case. 
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Fig.  8.  The  temporal  growth  rate,  as  calculated  from  the  dispersion  relation,  is  plotted  versus 
k  for  the  B0  —  5  kG  case. 
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Fig.  9.  Fourier  component  of  Bz(k  =  1.5  cm  *)  of  the  TE\\  mode  versus  time. 
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